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ABSTRACT 


To  support  SECEDE  II  data  analysis  and  interpretation,  EG&G 
(Bedford  Division)  is  analyzing  photographic  data  from  the  Spruce  Event. 
Fourier  analysis  of  digitized  data  frames  will  produce  power  spectra, 
autocorrelation  functions,  and  other  such  data  as  is  appropriate  for  event 
interpretation  and  correlation.  The  concept  of  Fourier  analysis  is  also 
being  examined  to  establish  the  validity  and  usefulness  of  tin  analysis 
results. 
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SECTION  1  , 

INTRODUCTION 

'  ;  ! 

.  ! 
i  i 

The  Bedford  Division  of  EG&G  hap  channeled  its  efforts  on  thd  SECEDE 

i  . 1 

program  in  two  directions.  First,  it  is  attempting  to  characterize  the  Struc¬ 
tured  environment  (i.  e. ,  striations)  created  in  a  late-time  barium  cloud.  ! 
Second,  it  is  acquiring  and  correlating  all  such  data  from  SECEDE  I,  II,  and 
III,  BIRDSEED,  certain  pre -SECEDE  releases,  and  fpom  other  environments. 
This  is  a  first  semi-annual  technical  report  describing  some  preliminary  work 
in  both  of  these  areas. 

•  *  1  I 

In  the  first  area,  EG&G  has  been  asked  to  proceed  with  a  digital  analysis 
of  appropriate  photography.  In  particular,  selected  frames  of  data  from 
Technology  International  Corporation  (TIC,)  photography  are  being  scanned  on 
a  Mann  Trichromatic  Microdensitometer,  and  the  density  and/or  radiance 

} 

profiles  generated  from  the  scanned  data  are  being  digitally  analyzed  using 
Fourier  techniques  to  produce  power  spectra,  auto -covariance,  and  auto¬ 
correlation  functions,  and  other  such  data  as  might  prove  useful,  i 

*  1  i 

In  parallel  with  the  analysis  of  photographic  data,  EG&G  is  also  under¬ 
taking  an  examination  of  the  applicability,  validity,  and  limitations  of  such' 
analysis.  R  has  become  apparent  that  members  of  the  SECEDE  community 
not  knowledgeable  in  the  areas  of  photographic  data  analysis,  digital  analysis,; 
and  Fourier  techniques  would  find  such  a  study  of  considerable  interest.  In 
particular,  EG&G  will  attempt  to  examine  ttye  realistic  and  practical  limits  of 
such  analysis  work,  and  will  attdmpt  to  give  the  user  a  basic  understanding  of 
the  meaning  of,  and  validity  of,  power  spectra,  Foufier  transforms,  and  other  ' 
such  data. 
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In  the  second  area,  EG&G  is  reviewing  all  available  SECEDE  docu¬ 
mentation  and  will  organize  pertinent  information  in  tabular  and/or  matrix 
formulations  rto  facilitate  correlation  and  interpretation  of  results.  It  is 
currently  anticipated  that!the  results  of  this  effort  will  be  presented  in  EG&G's 
final  report,  when  all  available  data  has  been  acquired. 

i  '  ' 

I 

Section  2  herein  presents  the  radiometric  data  which  is  being  Fourier 
analyzed.  Pertinent  frames  of  data  from  the  Spruce  Event  of  SECEDE  II  are 
identified  for  analysis  and  radiance  profiles  presented.  Section  3  begins  the 

derailed  examination  of  Fourier  analysis  as  applied  to  microdensitometer 
data. 

Appendix  A  contains  notes  prepared  during  the  SECEDE  Summer  Study. 
These  notes  detail  some  of  the  assumptions,  tools,  and  relationships  used  by 
EG&G  in  obtaining  radiometric,  information  from  photographic  data.  It  is 
reproduced  here  for  the  benefit  of  those  members  of  the  SECEDE  community 
who  did  not  receive  it  at  the  Summer  Study. 


f 


i 

i  i 
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SECTION  2 


SPRUCE  STRIA  TIONS 


The  EG&G  striation  analysis  effort  will  be  directed  at  Event  Spruce. 
In  particular,  four  frames  of  data  (Table  1)  will  be  analyzed  in  the  Fourier 
domain  for  spatial  frequency  content.  These  frames  were  selected  because 
they  best  bracketed  the  sphere  and  beacon  track  times,  and  because  the  data 

generated  could  potentially  be  used  to  correlate  with  the  sphere  and  beacon 
data. 


Table  1.  Event  Spruce  Striation  Analysis. 
T.  I.  C.  Film  No.  71723 


Frame 

Time  (after  release) 

146 

16  min,  5  sec 

164 

17  min,  52  sec 

185 

20  min,  28  sec 

193 

23  min,  36  sec 

The  referenced  frames  and  generated  radiance  profiles  (across  the 
striations)  are  shown  in  Figures  1-4.  One  hopes,  with  Fourier  analysis 
techniques  applied  to  the  digitized  radiance  profile,  that  the  small  scale 
structure  seen  superimposed  on  the  overall  cloud  background  can  be  charac¬ 
terized  in  a  power  spectrum  presentation.  These  profiles  are  currently  being 

digitally  processed;  the  results  of  the  analysis  will  be  presented  in  the  final 
report. 
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Figure  1.  Event  Spruce  as  seen  at  16  min,  5  sec  from  Tyndall. 
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Figure  2.  Event  Spruce  as  seen  at  17  min,  52  sec  from  Tyndall. 
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Event  Spruce  as  seen  at  20  min,  28  sec  from  Tyndall. 
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Figure  4.  Event  Spruce  as  seen  at  23  min,  36  sec  from  Tyndall 
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The  Spruce  ion  cloud  at  approximately  20  minutes  is  shown  in 
Figure  5  as  Frame  187;  from  T.I.C.  record  71721  (originally  a  color  record, 
printed,  here  ih  black-and-white).  This  record  was  exposed  from  Tyndall; 
and  assuming  an  ROA  (range  along  the  optical  axis)  of  300  km,  a  dimensional' 
scale  was1  placed  on  the  photograph  to  indicate  relative  cloud  and  striation  1 
dimensions.  ,  1  >  1 

1  i  ! 

Of  particular  concern  here  is  the  striated  Region  indicated  by  the 
arrow.  Are  these, characteristic  and/or  repetitive  striation  dimensions  and 
spacings,  which  can  be  measured  from  photography  such  as  this?  If  npt,  what 
data, can  be  gleaned  from  records  such  as  these,  and  what  does  this  data 
indicate  about  the  spatial  and  temporal  history  of  stria;tions  per  se.  Section  3 
begins  a  detailed  examination  of  this  question  and  presents  the  implications 
of  '(he  answers  to, that  question.  !  i 


i  i 

i. 
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Figure  5.  Event  Spruce  as  seen  from  Tyndall  on  T.  I.  C. 
record  71721  at  20  minutes. 


SECTION  3 


SOME  PRELIMINARY  THOUGHTS  ON  DATA  PROCESSING 

FOR  STRIA TIONS 


3.  1  INTRODUCTION 

Shortly  after  the  high- altitude  release  of  barium,  the  ion  cloud  is 
observed  to  form  striations  -  highly  structured  "tubes"  roughly  parallel  to  the 
geomagnetic  field.  The  striations  are  visual/ photographic.  They  map  the 
local  geomagnetic  field,  and  their  formation  is  the  result  of  interactions 

between  the  geomagnetic  field,  the  Afield,  ionospheric  winds,  and  the  ion 
cloud. 

In  order  to  interpret  the  available  SECEDE  data  and  to  predict  future 
data,  it  is  necessary  to  describe  the  striations,  at  least  statistically.  A 
najor  portion  of  the  SECEDE  effort  is  thus  directed  toward  developing  in¬ 
formation  on  the  growth  and  formation  of  striations. 

One  part  of  the  overall  striation  analysis  effort  is  directed  toward 
the  generation  of  spatial  frequency-domain  power  spectra  of  striated  ion-cloud 
radiance  profiles.  It  is  generally  believed  that  such  information  will  assist 
in  the  assessment  and  interpretation  of  SECEDE  data,  although  the  detailed 
connection  between  the  spatial  frequency  content  of  visible  striations  and 
the  various  signal  degradations  has  not  been  established.  Perhaps  the  con¬ 
nection  could  be  established  if  good  spatial  frequency  descriptions  were 
available. 

But  what  do  we  mean  by  "good"?  Ideally,  we  would  mean  descriptions 
permitting  an  exact  prediction  of  any  future  event;  that  is,  deterministic 
relations  or  equations  into  which  we  would  insert  the  various  experimental 
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parameters  (height  of  release,  quantity  of  free  barium,  etc.  ).  Those  equations 
would  then  be  solved  to  determine  such  characteristics  as  the  time  and  location 
at  which  striations  would  form,  their  size,  brightness,  and  distribution  in  space. 
Apparently  such  equations  are  presently  beyond  our  grasp.  In  fact,  the  SECEDE 
experiments  suggest  that  the  distribution  and  size  of  striations  are  the  result 
of  a  random  physical  process.  If  that  is  the  case,  completely  deterministic 
equations  are  unattainable. 

It  does  not  matter  whether  the  apparent  randomness  correctly  reflects 
the  fundamental  nature  of  the  process,  is  a  consequence  of  an  uncontrolled 
(random)  systematic  parameter  such  as  the  local  If-field,  or  simply  reflects 
our  lack  of  theoretical  understanding  of  a  deterministic  process.  Our  best 
current  theory  is  that  the  size,  location,  and  intensity  of  striations  are  the 
results  of  a  stochastic  process  for  a  given  event. 

What  we  can  mean  by  "good"  descriptions,  then,  are  those  from  which 
an  accurate  statistical  "prediction"  of  the  striations  which  might  be  observed 
in  some  future  event  can  be  made.  Suppose,  for  instance,  that  the  radiance 
profiles  of  individual  striations  have  the  same  waveform,  but  that  amplitude 
and  spacing  are  random.  In  that  case,  a  good  description  would  consist  of 
(1)  a  description  of  that  waveform,  (2)  an  amplitude  description  (at  least  the 
mean  amplitude),  and  (3)  spacing  statistics. 

3.  2  STATIONARITY  AND  ERGODICITY 

Stochastic  processes  are  analyzed  through  the  examination  of  various 
ensemble  averages;  that  is,  the  same  experiment  is  repeated,  in  principle, 
an  infinite  number  of  times,  and  the  results  of  each  are  treated  as  one  statistical 
sample  of  the  random  process.  However,  for  certain  restricted  classes  of 
stochastics  it  is  possible  to  obtain  the  desired  statistics  with  considerably 
less  than  an  infinite  number  of  identical  experiments,  perhaps  with  as  few 
as  one. 
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A  classification  tree  for  stochastic  processes  is  shown  in  Figure  6. 
Processes  classified  as  ergodic  have  the  property  that  the  statistics  of  a  single 
sample  are  identical  to  the  ensemble  statistics.  In  reality,  we  never  observe 
a  complete  sample  (since  it' s  infinitely  long)  but  rather  some  fraction  of  the 
complete  sample  which  might  be  called  the  sample  record.  From  one  such 
record,  we  can  obtain  an  estimate  of  the  statistical  parameters  of  the  process. 
The  accuracy  of  that  estimate  is  related  to  the  record  length.  That  accuracy 
can  be  improved  by  analyzing  a  longer  record  or  by  combining  data  from 
shorter  records. 


Figure  6.  Classification  of  Stochastic  Processes. 

3.  2.  1  Classification  of  Stochastic  Processes 

Figure  7*  shows  part  of  an  ensemble  of  sample  records.  If  we  read  the 
value  of  x  at  tj  for  each  record  and  divide  by  the  number  of  records,  we  obtain 
the  ensemble  average,  designated  mv( t. ),  at  t..  In  general,  »  (t. )  will  take  on 
different  values  as  tj  is  changed.  In  the  same  way,  we  may  define  the  ensemble 
average  autocorrelation  function,  R  (tj,  tJ#  +r),  which  will  generally  be  a 
function  of  both  tj  and  r.  The  equations  are 


*Many  illustrations  use  "t"  to  denote  the  independent  variable  and  x  as  the 
stochastic.  The  origin  of  most  of  this  work  is  in  communication  theory  where 
t  is  time.  The  reader  is  expected  to  make  the  appropriate  mental  substitutions, 
such  as  t-*x  and  x-»N(x). 
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If  these  averages  do,  in  fact,  vary  with  the  choice  of  tj,  then  the  process  is 
said  to  be  nonstationary.  But,  if  it  should  happen  that  the  ensemble  average 
and  ensemble  autocorrelation  are  independent  of  the  choice  of  t.,  the  process 
is  said  to  be  stationary.  *  In  that  case,  we  mav  write  u  (t. )  =  u  and 

It  is  also  possible  to  form  a  mean  and  autocorrelation  for  a  single  sample 
function: 


MxOO 


T 

tT.  ^  I  V1’ 


dt 


R x<r,  k)  « 


J. 

T  /  XjJt)  x^t  +  r)  dt 


where  k  denotes  the  k-th  sample  of  the  (stationary)  process.  If  both  Mx(k) 
and  Rx(r,  k)  have  the  same  value  for  all  k*  s,  then  the  process  is  said  to  be 
ergodic.  In  that  case,  the  various  sample  averages  and  ensemble  averages 
are  equal;  i.  e. ,  #1  (k)  ■  and  R  (r,  k)  ■  R  (r). 

XXX  X 


♦Strictly,  if  m  and  R  are  independent  of  t1#  the  process  is  ’’weakly  stationary". 

If  all  higher  order  moments  and  joints  moments  are  also  independent  of  tt,  the 
process  is  "strongly  stationary"  or  "strictly  stationary". 
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The  concept  of  stationarity  relates  to  the  ensemble  averaged, properties 
of  a  random  process.  In  practice,  however,  data  in  the  form  of  individual 
time-history  records  for  a  random  phenomenon  are  frequently  referred  to  as 
stationary  or  nonstationary.  A  slightly  different  concept  of  stationarity  is 
involved  here.  When  a  single  time-history  record  is  referred  to  qs  being 
stationary,  it  is  generally  meant  that  the  properties  computed  over  short  1 
time  intervals  do  not  vary  "significantly"  from  one  interval  to  the  nex{.  The 
word  "significantly"  is  used  here  to  mean  that  observed  variations. are  greater 
than  would  be  expected,  owing  to  normal  statistical  sampling  variations.  Hence 
the  single  sample  record  is  stationary  within  itself.  This  concept  of  sta¬ 
tionarity  is  sometimes  called  self- stationarity  to  avoid  yonfusion  with  the 

more  classical  definition.  i  , 

.  1  1  ! 

To  clarify  the  idea  of  self- stationarity,  consider  a  single  sample  record 

xfc(t)  obtained  from  the  k-th  sample  (function  of  random'  process  x(t).  Assume 
that  a  mean  value  and  an  autocorrelation  function  are  obtained  by  .time  averaging, 
over  a  short  interval  T  with  a  starting  time  oft.  as  follow's. 


t.i+T 


Mx(tl* 10 


T  f  Xk(t) 


dt 


tj+T 


Rx^l*4!  +T*k) 


4  / 


(t)  Xj^t  + 


t)  dt 


For  the  general  case  where  the  sample  properties  vary  significantly  as  the 
starting  time  t^  varies,  the  individual  sample  record  is  said  to  be  self- 
nonstationary.  For  the  special  case  where  the  sample  properties  do  pot  i 
vary  significantly  as  the  starting  time  t  varies,  the  sample  record  is  said  to 

A  i  , 

be  weakly  self- stationary.  If  this  requirement  is  met  for  all  higher  order 

» 

moments  and  joint  moments,  the  sample  record  is  s^id  to  be  strongly  self- 
stationary. 
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An  important  point  here  id  as  follows:  A  sample  record  obtained  from 
an  ergodic  random  process  will  be  splf- stationary.  Furthermore,  sample 
records  from  most  physically  interesting  nonstationary  random  processes 

I  !  '  : 

will  be  self-nonstationary.  Hence,  if  an  ergodic  assumption  is  justified 
(as  it  is  for  most  adtual  stationary  physical  phenomena),  verification  of 
self  -  stationarity  for  a  single  sample  record  will  effectively  justify  an  assump¬ 
tion  of  stationarity  ^ind  ergodicity  for  the  random  process  from  which  the 
sample  record  is  obtained. 

1  :  l 

Ergodic  random  processes  are  clearly  an  important  class  of  random 
processes  since  all  properties  of  ergodic  random  processes  can  be  deter¬ 
mined  by  perfbrming  time  averages  over  a  single  sample  function.  Fortunately, 
in  actual  practice,  random  data  representing  stationary  physical  phenomena 
are  generally  ergodic.  It  is  fpr  this  reason  that  the  properties  of  stationary 
random  phenomena  can  be  measured  properly,  in  many  cases,  from  a  single 

,i  1 

observed  time -hi story  record. 

3i.  2.  2  Nonstatiqnary  Stochastics 

.  t  i* 

Nonstationary  random  processes  include  all  random  processes  which 
do  not  meet  the  requirements  for  stationarity.  Unless  further  restrictions 

■  i 

are  imposed,  the  properties  of  nonstationary  random  processes  are  generally 

time-varying  functions  which  can  be  determined  only  by  performing  instantaneous 

‘  «  *  * 

averages  over  the  ensemble  of  sample  functions  forming  the  process.  In 
practice,  it  is  often  not  feasible  to  obtain  a  sufficient  number  of  sample  records 

i  I 

1 

to  permit  the  accurate  measurement  of  properties  by  ensemble  averaging. 

This  fact  has  tended  to  impede  the  development  of  practical  techniques  for 

)  1 . 

measuring  and  analyzing  nonstationary  random  data. 

In  fact,  a  totally  adequate  methodology  does  not  as  yet  exist  for  the 
analysis  of  all  types  of  nonitationary  data,  partly  because  of  the  fact  that  a 
nonstationary  conclusipn  is  generally  a  negative  statement  specifying  the 
lack  of  stationary  properties,  rather  than  defining  the  precise  nature  of  the 

I  I  ! 
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nonstationarity.  On  the  other  hand,  when  a  process  is  deemed  stationary, 
certain  positive  results  are  known  which  apply  to  all  stationary  data.  For 
nonstationary  data,  special  techniques  must  be  developed  which  apply  only 
to  limited  classes  of  these  data. 

Illustrated  in  Figure  8  are  the  three  basic  and  most  important  categories 
of  nonstationary  data:  (a)  time-varying  mean  value,  (b)  time-varying  mean 
square  value,  and  (c)  a  combination  of  (a)  and  (b).  As  we  shall  see,  the 
radiance  profile  of  the  striated  ion  cloud  is  probably  nonstationary  (category  c); 
however,  it  may  be  possible  to  treat  the  striations,  per  se,  as  ergodic. 


Tim# 


Tim# 


Tim# 


Figure  8 .  Examples  of  nonstationary  data. 


3.3  SECEDE- TYPE  STRIA  HON  DATA 

Figures  1  through  4  contain  typical  radiance  profiles  for  striated 
barium  ion  clouds.  To  the  extent  that  these  scans  are,  in  fact,  typical,  it 
is  evident  that  the  ion  cloud  radiance  profile  data  has  a  distance  varying  mean 
and  mean  square.  (Compare  with  Figure  8. )  The  data  is  nonstationary  (thus, 
nonergodic).  We  are  not,  however,  particularly  interested  in  the  ~r^s 
behavior  of  the  ion  cloud,  nor  do  we  care  about  gradation  in  the  general  sky 
brightness.  If  we  could  subtract  these  two  effects*  from  the  radiance  profile, 
we  might  be  left  with  a  data  record  which  is  ergodic. 

To  the  best  of  our  knowledge,  such  a  subtraction  has  not  been  attempted 
for  SECEDE-type  striation  data;  thus,  we  do  not  know  exactly  how  to  proceed. 
We  are  inclined  to  try  simple  procedures  first,  and  the  remainder  of  this 
section  will  be  concerned  with  one  fairly  simple  approach  which  also  has 
considerable  flexibility;  that  is,  to  subtract  a  running  average  from  the  data 
record.  The  average  should  be  centered  on  the  data  point  being  operated  upon. 
This  process  is  essentially  high-pass  filtering,  with  a  "turn-on"  frequency 
determined  by  the  averaging  interval.  The  data  of  Figures  1  through  4  suggests 
that  the  background  radiance  level  (basically  interpreted  to  include  gross  ion- 
cloud  radiance,  sky  brightness,  and  striation  DC)  occurs  at  much  lower  fre¬ 
quencies  than  the  striations,  so  this  simple  subtraction  process  may  provide 
satisfactory  results. 


It  will  also  be  necessary  to  suppress  the  DC  term  required  for  physically 
meaningful  radiance.  That  suppression  is  desirable  for  other  reasons  dis¬ 
cussed  later.  If  desired,  the  DC  power  can  be  added  directly  to  the  spatial 
frequency  power  spectra. 
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A  Cautionary  Note 


A  literal  high- pas  ‘  filter  applied  to  the  Fourier  transform 
will  not  give  equivalent  results  to  those  obtained  using  the 
subtraction  process  We  are  working  with  sampled, 
digitized  records  of  finite  length.  In  consequence,  the 
numerical  Fourier  transform  at  each  frequency  contains 
a  contribution  from  all  other  frequencies.  The  numerical 
transform  is  the  convolution  of  ideal  transform  and  a  sinc- 
function*  corresponding  to  the  length  of  the  data  set.  The 
sine -function  propagated  by  DC  will  affect  all  higher  fre¬ 
quencies  and  may,  in  fact,  dominate  the  power  spectrum 
(as  is  discussed  further  herein). 

The  effect  of  the  convolution  occurs  simultaneously  with 
the  transform.  It  cannot  be  unfolded.  If  high  frequency 
data  has  been  modified  substantially,  a  high-pass  filter 
cannot  retrieve  the  ideal  data. 

Basically,  what  we  are  doing  is  writing  the  total  signal  N(x)  as  the  sum 
N(x)  *  B(x)  +  S(x) 

of  a  nonstationary  background,  B(x),  which  is  of  little  or  no  interest,  and  a 
dc  suppressed  striation  term,  S(x),  which  we  hope  to  treat  as  ergodic.  For 
a  single  record,  the  Fourier  transform  is 

N(w)  *  S(w)+'S(w) 

so  that  the  ideal  power  spectrum  is 

P(N,  w)  =  N(  2)  N(u)  *  P(B,  w)  +  P(S,  u>)  + 

[cross  terms  in  B  and  S] 


* 


sine  (x)  =  (sinx)/x 
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If  B  and  S  are  disjoint  in  frequency  space,  the  cross  terms  are  identically 
zero.  The  ideal  power  spectrum  of  a  single  record  would  simply  be  the  sum  of 
the  power  spectra  of  B  and  S.  Since,  however,  the  actual  spectrum  contains  a 
convolved  sine -function,  it  is  not  possible  to  simply  separate  P(S,  <4  out. 

On  the  other  hand,  a  prescription  can  be  given  which  hopefully  separates 
B  from  S  directly  and  which  essentially  guarantees  that  the  resulting  functions 
are  frequency  disjoint.  Since  P(B,  <4  and  P(S,  <4  can  be  formed  from  the 
separated  functions,  it  appears  that  no  information  on  the  total  power  spectrum 
of  the  record  is  lost.  However,  one  is  able  to  treat  the  interesting  part, 

P(S,  w),  separately. 

3.  3.  1  Ergodicity  of  S(x) 

Suppose,  then,  that  the  separation  is  effected.  One  must  ask  whether 
anything  has  been  accomplished;  that  is,  does  the  k-th  record  power  spectrum 
Pk@>'  <4  tell  us  anything  about  the  power  spectrum  of  striations  as  a  general 
class.  Once  again,  the  answer  isn't  clear.  It  depends  on  whether  some 
apparently  reasonable  assumptions  are  correct. 

If  we  assume  that  the  radiance  profiles  of  individual  striations  have  the 
same  form  and  that  their  location  and  amplitude  are  random,  then  it  can  be 
shown  that  the  striation  record  waveform,  S,  is  ergodic.  And,  in  fact,  if  v 
is  the  striation  ensemble  average  waveform,  then 

P(S,  <4  oc  P(v,  <4 

We  must  be  careful  in  interpreting  this  expression  since  we  are  dealing  with 
statistical  quantities.  Once  the  constant  of  proportionally,  k,  has  been 
determined,  P(S)  and  kP(v)  are  statistically  indistinguishable.  But  that  does 
not  mean  that  they  are  numerically  identical  at  each  wavelength,  u>.  As  the 
length  of  the  record  s(x)  increases,  the  statistics  improve,  which  is  to  say  the 
variance  [(P(S)  -  kP(v))2]  decreases.  Only  in  the  limit  of  an  infinite  record 
length  is  the  variance  of  power  estimates  at  each  wavelength  reduced  to  zero; 
thus,  only  in  that  limit  is  a  numerical  identity  obtained.  For  records  of 
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finite  length,  the  variance  of  the  estimate  at  each  wavelength  can  be  determined. 
Thus,  the  power  spectrum  of  an  ergodic  record  of  finite  length  gives  a,n  estimate, 
with  known  variance,  of  the  ensemble  average  power  spectrum. 

*  '  I 

To  illustrate  how  the  constraint  of  proportionality  is  obtained,  cpnsider,  a 
restricted  case  of  the  above  assumption  in  which  the  individual  striation  Radiance 
profile  amplitudes  are  fixed  and  identical.  The  stHations,  however,  occur 

*  •  I 

(replicate)  with  random  spacing.  Let  the  waveform  to  be  replicated  be1  A(t). 

(In  sketching  the  waveforms^  A(t).will  be  taken'to  be  a  square  pulse.  The 
formalism  is  general,  however.  )  The  Fourier  transform  of  the  replication 
element  is  A(u>)  (Figure  9)  •  ,  , 


A(u>)  .= 


,7S  f* 

T- eo  * 


A(t),  e  ^  dt 


and  is  in  general  a  complex  quantity.  The  power  spectrum  of  the  waveform  is 


(w) !  =  AM  A(u) 


which  is  a  real  function  of  w.  The  phase  relations  between  the  variojus  wave¬ 
lengths  has  been  lost. 


A(t) 


n 


-»>  t 


Figure  9.  A  single  feature  and  its  trarisform1.  ' 

'  i 

I 

Now,  let  a  replication  bf  A(t)  occur  at  a  later  time.  The’ waveform  is 
B(t).  The  total  waveform  is  A(t)  +  B(t),  ,as  illustrated  in  Figure  10.  Since 
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j  '  I  1  1  ,  1  1  I  . 

Fourier  transformation  is  a  linear  process,  the  transform  of  the  sum  is  1 

.  < 


; 


Thus,  the  transforms  of  B  and  A  are  identical  apart  from  phase  factors. 
Clearly  then,  if  the  total  waveform,  Z(t),  consists  of  replication  of  A (t)  at 
tines  characterized  by  tj,  t2,  ...  t  (see  Figure  11). 


Z(L) 


ia  nn 

n 

s - - - 

- »• - — » 

— 1 — i 

it 


Figure  11.  The  total  record  of  randomly  replicated  features. 
The  transform  is 

ZM  .  AM  [1  +e'i‘<l  +  e-,v*2  +  . . .  +  p-Hi  ] 
ami  the  power  spectrum  of  Z(t) 

Z*Z  .  A«A  [  1  +  e'il4>  +  e-,“,2  +  . . .  +  ] 


[1  +  c1"*1  +  etu*2  +  . . .  +el“»n] 


vhore  we  have  defined  tQ  •  0  so  that  e*iwt°  •  1.  Since  the  sums  are  over 
different  indices 


[N  N 
j-o  i«0 


Z*Z  ■  A*A 


23. 


We  can  conaider  that  the  aura  la  over  a  two- dime nalonal  array  of  the 
( tj.  tj } .  where,  for  inatance,  "i"  la  a  column  index  and  "j"  a  row  Index. 


Break  the  double  aum  into  three  plecea,  one  piece  conaietlnf  of  all  terma 

on  the  diagonal,  the  aecond  all  terma  above  the  diagonal,  and  the  third  all 
terma  below  the  diagonal 


2*2  *  A*A 


[r 

L  i.j.o 


-lu(t. -t,)  "  -lult.-t.) 

•  +2-  •  i 


fj.-vvl 

<  J 


Each  element  of  the  fir  at  ^  ia  identically  I,  and  the  aecond  and  third 
terma  are  complex  coojugatea.  Thua 


~  ~ -  -iu(t,  - 1.) 

2*2  ■  A*A  N+  2 Re  ^  e  J  1 

J>i 
J  *  0 


m'm 

2*2  *  A*A  N  ♦  2  £  coew(tj  -tj] 
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To  this  point,  the  derivation  is  rigorous.  If  the  time  at  which  the 
replication!  occur  ia  random,  the  argument  of  the  coeine  is  alao  random. 

If  N  ia  large  enough,  one  expects  to  aee  as  many  negative  terma  aa  poaitive 
terma  in  the  remaining  summation.  Thus,  if  N  is  large  enough,  the  remaining 
summation  goes  to  aero  and 


z.z  -  n;.; 

N large 

In  words,  if  enough  randomly  spaced  replications  are  included,  the 
power  spectrum  of  the  total  waveform  ia  simply  N-timea  the  spectrum  of  a 
replication  unit. 


To  complete  the  analysis,  consider  that  the  real  waveform  consists  of 
randomly  timed  replications  of  A(t)  plus  some  other  signal  C(t): 


Then 


N 

Z(t)-C(t)+  5^A(t-tt> 
i»0 


and 


~  -itft.  ^  *  -Hut. 

Z*Z  ■  C*C  +  NA*A  ♦  C*A  2J  •  +CA*  e  J 

i*0  j-0 

The  last  two  terms  are  complex  conjugates,  and  their  total  contribution  is 
just  twice  the  real  part  of  either  terra.  In  this  case,  it  is  not  possible  to 
write  the  real  part  in  a  simple  form;  however,  it  is  possible  to  write  a 


superficially  simpler  form  as  a  basis  to  argue 
Let 

D  ■  C*  A  *  d  ♦  i6  and  Q»^ 

i 

then  the  last  two  terms  contribute 

2  Re  D  G  ■  2  (df  +  6  a) 

since 

«■  «r0 

and 

a  •  ^sin  wtj  »  0 
By  the  same  argument,  the  last  terms  are  negligible  to  first  order  and 

<*«  "<  «w 

Z*Z  -  C*C  +  NA*A 

In  words,  the  power  spectrum  of  a  waveform  containing,  among  other 
things,  N  randomly  spaced  replications  of  a  particular  signal  (A(t))  will 
contain  (to  first  order)  N  -times  the  power  spectrum  of  A. 

The  more  general  assumption  which  we  wish  to  make  about  striations 
is  that  the  amplitude  is  also  a  random  variable.  (This  is  sometimes  termed 
Impulse  Noise". )  In  that  case,  S(x)  is  the  sum  of  many  similarly  shaped 
transient  pulses  (see  Figure  12). 


that  the  term  is  negligible. 


-iut. 


8  + 1*  . 
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Figure  12.  Sample  functions  x  (t)  for  five  examples  of  random  processes. 


26- 


00 

s<*>  *  2  ^  v(x-*k) 
k»l 

whole  ihape  ii  given  by  v  v  (x).  with 
00 


v(x)  e"iwx  dt  ■  V  (w) 


/• 

*'-00 


while  the  pulse  amplitude  afc  is  a  random  variable  with  finite  variance, 
and  the  x's  are  random  locations  determined  by  the  state  changes  of  a 
Poisson  process  with  mean  count  rate  e.  The  process  is  stationary 
and  ergodic  if  started  at  x  ■  -  one  has 

00 

E{S(x}|  ■  »E  (ak)  I  v(x)  dt 

-00 

00 

Rg(x)-M+  «E  (ak)  I  Vix)  Wx  +  1)  dt 

.00 


P(S,«)  *  (u)  ♦  jL.  e  {^  )  |  V(w)| 


where  E  (u  )  denotes  the  expected  value  of  u. 


APPENDIX  A 


RADIOMETRIC  INTERPRETATION 
OF 

PHOTOGRAPHIC  RECORDS 


A.  1  INTRODUCTION 

These  notes  present  a  brief  outline  of  the  mathematical  formalism 
and  the  concepts  used  by  EG&G  to  obtain  source  radiance  and  power  infor¬ 
mation  from  photographic  records.  The  notes  were  prepared  during  the 
1971  SECEDE  Summer  Study. 

A.  2  BASIC  CONCEPT 

Radiometric  interpretation  of  photographic  records  is  based  on  the 
assumption  that  the  measured  density  of  an  element  of  film  is  a  monotonic 
function  of  the  effective  exposure,  as  defined  by  Eqs.  (1)  and  (2),  provided 
that  other  variables  which  affect  density  are  held  constant  (e.  g. ,  exposure 
duration,  processing,  etc. ).  Therefore,  if  a  calibration  curve  (a  D-log  u 
curve)  is  generated  by  measuring  the  densities,  Dj,  of  known  exposures, 

1 1*  on  film  of  the  same  emulsion  as  (and  processed  with)  the  data  film, 
the  effective  exposure  on  the  data  film.  Eg  (x,  y),  associated  with  a 
measured  density,  D{  (x,  y),  may  be  obtained  from  this  curve. 

A.  3  DENSITY  DETERMINATION 

Density  measurements  may  be  made  with  any  of  the  sever*,  com¬ 
mercially  available  microdensitometers,  with  the  understanding  that 
calibration  densities,  Dj,  and  image  densities,  D{  (x,y),  be  measured  with 
the  same  microdensitometer,  using  the  same  optics,  and  as  close  together 
in  time  as  is  feasible.  This  will  eliminate  all  problems  created  by  specular 
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versus  diffuse  measurements*  and  reduces  the  measurement  to  a  comparison 
one  in  which  the  calibration  and  image  densities  are  used  only  as  transfer 
measurements  of  exposure.  •  !  ( 

!  : 

A.  4  CALIBRATION  PROCEDURE 

i  •  i 

A.  4.1  Step  Wedge  Exposure 

.  i  •  i  i 

As  a  minimum*  each  data  record  should  receive  four  step-wedge 
exposures  in  a  calibrated  sensitometer  (two  pre-shot  and  two  post-shot), with  ■ 
each  pair  placed  head-to-head.  Sensitometer  exposure  times  should  be  chosen 
to  most  nearly  correspond  to  data  film  exposure  times.  The  step  wedge 
exposures  should  be  processed  at  the  same  time  as,  and  adjacent  to,  the 
associated  film  record. 

A.  4. 2  D-Log  E  Curve  , 

The  effective  exposure,  E^,  for  the  i^  step  of  the  step  wedge  is 
taken  to  be  1  , 

/ 

-  J  u(A)  Tt(A)  S(\>  d(\ )  ergs  ,  '  (1) 

0  cm^  i 

i 

»  : 

i  .  i 

where 

‘  '  i 

u(A)  5  spectral  energy  depsity  on  the  sensitofoeter 
platen  (ergs/cm^-A), 

T  (A)  *  dimensionless  diffuse  transmittance  of  the  1 

ith  step,  and  , 

I 

S(X)  3  dimensionless  film  sensitivity.  ' 

EG&G  generally  scans  each  of  the  four  step  wedges  three  times  along 
parallel  paths  to  produce  twelve  sets  of  density  data  (thirty-six  for  color 
film).  Trial  D-log  E  curves  are  constructed  for  each  of  the  twelve  scans, 
and  these  are  reviewed  in  order  to  make  adjustments  and  to  eliminate 

1  '  !  . 

I 

I 

I  > 
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inconsistent  data.  A  final  b-log  ^  curve  is  generated  as  an  arithmetic 
qiean  of  the  acceptable  trail  curves. 

i 

I 

A.  5  RADIANCE,  'POWER,  AND  ENERGY  EQUATIONS 
A.  5. 1  Radiance 

,  —  |  i 

i  For  an  imaging  optical  system  focused  at  infinity,  the  effective 
exposure  at,  (x,  y)  on  the  film  is  given  by* 

7  . 

,  EI  (x,y)  *  l|g£  '  f  N(X)S(X)T(X)dX  ergs,  (2) 

4^7^  0  cm’ 

f  • 

i  , 

where  1 

i  i 

!  N<x)  =  spectral  radiance  of  point  (X,  Y)  in  the  object 

'  plane  corresponding  to  point  (x.y)  in  the  image 

plane  (watts/cm2<-ster-A), 

'  i  , 

T(X)  *  the  total  spectral  transmittance  along  the  optical 
1  path,  i  1 

(f/)  =  the  aperture  of  the  optical  system  (focal  length 
divided  by  the  diameter  of  the  lens  aperture),  and 

ieXp  =  the  exposure  time  at  point  (x,y)  (seconds). 

To  unfold  Eqn.  (2)  for  N(X)»  assume  that  the  relative  shape,  N(X). 
of  N(X)  is  known  but  not  its  magnitude,  NQ.  Then, 

t  1  N(X>,  r  Nq  N(X),  ,  (3) 


and  Eqn.  (2)  yields 


N(X)  'HhmL  1 

x  t 

i  exp  i 


A 

N(X)  _  .  . 

-  E! 

J  N(X)  T(X)  S(X)  dX 


Eqn.  (2);  assumes  that  the  source  area  considered  is  on  or  near  the 
'Optical  axis. 


i 


i 


i 


The  total  radiance  in  the  wavelength  range  (a,  b)  is  then 


Nab<X.Y>.  J  N(X)  dX 


r 

=  10  ^  4 (f /V^  a*^  N(X)  dX  Ejix,  y)  watts 

*  *exp  f  N(X)  TOO  S(X)  dX  cm^-ster 


Note  that  TOO  may  be  written  as 


T(X)  =  Ta(X)  Tw<X)  Tp(X)  Tp(X)  Tl(X). 


where 

TaO)  =  spectral  transmittance  of  the  atmosphere, 

TW<X)  spectral  window  transmittance, 

Tp(X)  =  spectral  filter  transmittance, 

Tp(X)  =  spectral  prism  transmittance,  and 
Tl(X)  =  spectral  lens  transmittance. 

For  black  and  white  film,  the  limits  of  integration  (a.b)  are  usually  3800  - 
6800  A.  For  color  film,  each  layer  is  treated  separately  with  the  limits 
usually  taken  as  3800  -  4800  A,  4300  -  5800  A,  and  5800  -  6800  A. 

A.  5.  2  Power 

Assume  that  the  radiating  source  is  located  on  the  optical  axis  and 
that  its  dimensions  are  small  compared  to  the  distance  between  the  source 
and  the  camera.  Then,  if  the  source  is  (1)  an  optically  thin  volume  rather 
or  (2)  a  spherical  Lambert  surface  radiator,  the  total  radiated  power  in  the 
wavelength  interval  (a,  b)  is  given  by 


P  .  =  4ir 
ab 


J  j  Nab  <X'Y> 


dX  dY 
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or, 


P  =  4ir 
ab 


(tft)  J  J\ b<x>y,dxdy 


watts 


(8) 


where  the  integration  is  performed  over  a  designated  area  of  the  film 
(perhaps  the  entire  frame),  and 

ROA  =  range  along  the  optical  axis  and 
FL  =  lens  focal  length. 

Once  Nafa  (X,  Y)  has  been  evaluated  from  Eqn.  (5),  may  be  calculated 
from  Eqn.  (8).  When  the  above  assumptions  do  not  hold,  special  calculation 
procedures  must  be  followed. 

A.  5.3  Total  Energy 

Through  a  consideration  of  successive  frames  of  a  record,  the  time 

dependence  of  Pab  may  be  determined.  Thus,  the  total  radiated  energy  in 

the  range  (a,  b)  and  in  the  time  interval  (t0-t  )  may  be  written  as 

t  2  1 

2 

Pab  ^  dt  30ules-  (9) 

Equations  (5),  (8),  and  (9)  are  the  basic  equations  employed  in  obtaining 
radiometric  information  from  photographic  records. 

A.  6  FACTORS  IN  THE  EQUATIONS 

A.  6. 1  Exposure  Time 


The  exposure  time  dexp)  should  be  measured  directly  or  computed 
from  measurements  of  film  velocity  and  camera  characteristics. 

A.  6.  2  Lens  Spectral  Transmittance /Optical  Aperture 


EG&G  obtains  the  lens  spectral  transmittance  (T  (X))  from 

Li 

measurements  of  the  actual  lens  at  the  nominal  setting.  The  curve  thus 
generated  is  not  the  true  spectral  transmittance  of  the  lens,  but  includes 
any  differences  between  the  nominal  f/  and  the  effective  f/.  The  curve 
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generated  is 

T  (X)  =  T,(X)  {tl)l  n0mtn*1  . 

(f/)  effective 

where  T^X)  is  the  true  spectral  transmittance  of  the  lens.  The  nominal  f/ 
is  used  in  the  radiance  equation. 

A.  6. 3  Atmospheric  Spectral  Transmittance 

The  atmospheric  spectral  transmittance  (T^(X))  can  be  a  unique 
quantity  for  each  photograph.  Ideally,  measurements  of  TA(X)  should  be 
made  for  each  exposure,  or,  since  this  is  almost  an  impossible  task, 
a  mathematical  model  of  atmospheric  transmission  can  be  constructed  so 
that  Ta(X)  can  be  calculated  as  both  a  spectrally  and  spatially  varying 
quantity. 

A.  6.  4  Spe  ^al  Trans mittances 

The  several  spectral  trans mittances  used  in  Eqns.  (5)  and  (6) 
should  all  be  measured  values. 
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